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ABSTRACT 


We present an approxiniation framework for computation (in finite dimensional 
spaces) of Riccati operators that can be guaranteed to converge to the Riccati 
operator in feedback controls for abstract evolution systems in a Hilbert space. 
It is shown how these results may be used in the linear optimal regulator problem 
for a large class of parabolic systems. 
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1. Introduction 


In this note we consider feedback controls for parabolic partial differential 
equations and the related Riccati operator theory when an infinite horizon 
integral quadratic cost functional is optimized. A general convergence frame- 
work for approximation ideas which can be used in computational techniques is 
developed in the context of the regulator problem theory pursued by Gibson in 
several recent investigations [ 8, 9, 10]. To illustrate our ideas we shall 
consider a specific model problem: The infinite horizon regulator problem for 

the parabolic control system 



for t > 0, X € C r'^, with Dirichlet boundary conditions 0 and known 

initial data yl^^Q = Consideration of this model problem is motivated by our 

desire to develop efficient computational schemes for optimal control problems 
in connection with the insect dispersal investigations detailed in [1, 2]. 

These problems (involving parabolic pav'tial differential equations) will entail 
distributed controls (e.g., spraying of pesticides over a region with frequency 
and intensity of spraying constituting important control variables). We expect 
the theoretical results presented in this paper to form a sound foundation for 
the development in the near future of computational procedures for feedback 
controls in such problems. 

In section 2, we state carefully a convergence theory for approximate 
Riccati operators that is essentially a modification and refinement of the theory 
presented by Gibson in [8]. (In an appendix, we indicate details as to how our 
framework follows from the results of Gibson.) We then in section 3 state 
precisely our control problem for the system (1.1) and show that under reasonable 



2 


assumptions (which imply a certain "preservation of exponential stability under 
approximation" condition) the abstract framework of section 2 can be used to 
guarantee convergence of approximate solutions in the event the basic approxi- 
mation scheme enjoys rather fundamental convergence properties. These are 
sufficiently relaxed to allow a generous number of practical schemes (modal, 
splines of several types and orders) to fall within our treatment. 

A concluding section contains remarks on the potential usefulness of the 
results in this presentation. 
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2_^ Approxin iati on of an Abs tra ct Linear Regulator Problem 

In this section we summarize approximation results for an abstract linear 
optimal regulator problem that we shall subsequently employ in our treatment 
of parabolic systems. The results given here involve a minor but important 
modification of the abstract theory developed by Gibson in [8], Specifically 
our presentation is formulated so as to facilitate approximation of the regulator 
problem by a sequence of finite dimensional state space problems, each defined 
on a subspace of the state space of the original problem. Gibson's presentation 
[8] requires the approximating problems each be defined on the entire original 
state space and as we shall explain below, this can lead to some tedious technical 
considerations. Our modified framework of this section really follows directly 
from that of Gibson, but we shall defer to an appendix a detailed explanation 
of this aspect of our considerations. 

We suppose throughout that H and U are Hilbert spaces, that A : domAC H H 
is the infinitesimal generator of a strongly continuous or Cq semigroup T(t) 
on H and that Be C(U,H). We consider a control system in H given by 

y(t) = Ay(t) +• Bu(t), t > 0 , 

(2.1) 

y(0) = yQ , 

and an associated performance measure 

a> 

(2.2) J(yn,u) = / {<Dy(t) ,y(t)> + <Qu(t) ,u(t)> }dt 

^ 0 

where D € C(H), Q€ c(U) are selfadjoint and satisfy D > 0, Q > 0. Our 
fundamental abstract linear optimal regulator problem can then be stated as 

2 

(R) Minimize J(yQ,u) over u e L (0,'";U) subject to y = y(*;u) 
satisfying (2.1). 
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2 

We shall say that a function uG L (0,«;U) is an admissible control for the 
initial state Yq € H if J(Yq,u) is finite. As usual, a certain algebraic 
Riccati equation will play a fundamental role in our analysis and an operator 
n € C(H) is called a solution of the algebraic Riccati equation (A.R.E) if n 
maps domAinto doniA*and satisfies on H the equation 

(2.3) A*n + iiA - + D = 0. 

Here A* is the Hilbert space adjoint of A and we recall [ 4, p. 51] that it is 
the generator of the Cq semigroup T(t)* which is adjoint to T(t). We note that 
if IT satisfies (2.3) on domAthen (2.3) can be taken as an equation on H since 
nBQ ^B*]i - D is bounded so that A*ii + [[A has a bounded extension to all of H. 

The following result is taken from [ 8, Theorems 4.11, 4.6]. 

Theorem 2. 1 . Let A, B, Q, D be as given above. Then there exists a nonnegative 
selfadjoint solution n of the algebraic Riccati equation (2.3) if and only if 
for each yg e H, there exists an admissible control. If this latter holds, 
then the unique optimal control and corresponding trajectory for (R) are given 

by 

(2.4) u(t) = - Q"^B*nJ(t) 

(2.5) y(t) = S(t)yg , 

where is the minimal nonnegative selfadjoint solution of the A.R.E. (2.3) and 
S(t) is the Cg semigroup generated by A-BQ~^B*n^. If |y(t;u)| 0 as t “ for 

any admissible control (this is guaranteed for example by the condition D > 0), 
then is the unjque nonnegative selfadjoint solution of the A.R.E. If 
D > 0, then we also have that S(t) is uniformly exponentially stable. 
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In this theorem the term minimal for a selfadjoint operator is in reference 

to the usual ordering of selfadjoint nonnegative operators on a Hilbert space. 

We note that the iiiinimal solution of (2.3) can be obtained as the limit of 

a sequence of Riccati operators for associated finite interval regulator problems 

(see [ 8, Theorems 4.10, 4.11]) in a manner analogous to the usual procedur-e 

for finite dimensional state space regulator problems [13]. 

We next formulate a sequence of approximate regulator problems and present 

N 

a convergence result for the corresponding Riccati operators. Let H , 

N = 1,2,..., be a sequence of finite dimensional closed linear subspaces of 

H and : H >■ be the canonical orthogonal projections. Assume that T^(t) 

N N / Nx 

IS a sequence of Cq semigroups on H with infinitesimal generators A (E C(H ). 

Given operators C. l’(U,H^) and D^€ C(H*'^) we consider the family of regulator 

probl ems: 

(R^) Minimize J^(y^(0),u) over uG L^(o,'^>^';U) , 


where 


( 2 . 6 ) 


y^t) = AV(t) + B%(t), 

y^O) = y^ = , 


t > 0 , 


and 

(2.7) j'^(y^(0),u) = / {<DV(t),y^(t)> + <Qu(t),u(t)>}dt. 

0 

We note that since : U ->■ the trajectories of (2.6) evolve in and 

consequently (R^) is a linear regulator problem in the finite dimensional state 
N 

space H so that finite dimensional control theory is applicable here. We 
shall need several assumptions in a convergence statement regarding solutions 
of (r''^) and (R). 
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(HI): For each Yq G there exists an admissible control u^C L^(0,«>;U) for 

(R^) and any admissible control for (2.6), (2.7) drives the state of 
(2,6) to zero asymptotically. 

(H2): (i) For each z G H, we have T*'*(t)P*'^z T(t)z with the convergence 

uniform in t on bounded subsets of [0,«). 

(ii) For each z G H, we have T*^(t)*P^z -> T(t)*z with the convergence 
uniform in t on bounded subsets of [0,«>). 

(iii) For each v G U, B^v > Bv and for each z € H, B*'^*z -> B*z. 

(iy) For each z C H, D^p'^z > Dz. 


We remark that (H2)(i) implies in particular (take t = 0) that P*^z -> z 
for each z 6 H and in this sense we have the subspaces approximate H. 

If assumption (HI) holds, then the optimal control for (R^) is given 
in feedback form by 

(2.8) u^(t) = - Q“^B'^*iiV(t) 


where G C(H^) is the unique nonnegative selfadjoint solution of the 


algebraic Riccati equation on H 


N 


(2.9) + II V - = 0, 


and y^ is the corresponding solution of (2,6) with u = u^. Moreover 

,N, N -Ny „ „N N ,N 
J (yg>^ ) * 


We also have the following fundamental convergence results. 
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Theorem 2.2 . Suppose (HI), (H2) hold, Q > 0, D > 0 and ^ 0 and let li*'* denote 

N 

the unique nonnegative selfadjoint Riccati operators on H for the problems 
N\ 

(E ). Further assume that a unique nonnegative selfadjoint Riccati operator 
on H for the problem (r) exists. Let S(t) and be the semigroups generated 

by A - BQ'^B*n and - B^Q"'b^*IT^ on H and H^, respectively. If there are 

positive constants , Mg and u independent of N and t such that 

(2.10) |s'^(t)l for t > 0, N = 1,2,..., 

' 

and 

(^-n) l^N - ^^2 ’ 

then 

(2.12) > iiz for every z € H, 

(2.13) S^^(t)P^'*z >- S(t)z for every z e H, 

where the convergence is uniform in t on bounded subsets of [0,'^>)> 
and 

(2.14) |S(t)| < M^e’"’'*^ for t > 0. 

We present a proof of Theorem 2.2 in the Appendix. Meanwhile we remark 

N N N / Nx 

that under the hypotheses of this theorem, n P is an extension of n e C(H ) 

to all of H. If D^, A^^ are replaced by A^V^, respectively and (2.9) is 

N N 

considered as an equation on H, then n P is its unique minimal nonnegative 
selfadjoint solution. 

Theorem 2.2 is essentially contained in [8]. The main difference between 
the theorem here and the result in [8] is, as stated earlier, that here each 
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of the finite dimensional state problems (R^) is defined in the subspace 
only, whereas in [8], Gibson requires that the approximate regulator problems 
be defined on the entire space H. This causes some unnecessary technical 
difficulties: First note that if D > 0 and = P*'^D (as an operator in H^), 

then > 0 on But = 0 on This difficulty can be circumvented by 

considering instead = P^D + I - P^ as an operator on H -- see [ 9, p. 698]. 

To explain a second disadvantage to the formulation of (R^) on all of H, 
let us assume that |T^(t))^l^ _< Me~“^ for positive constants M and a. This 
allows one to infer existence of Riccati operators on however, if the 
semigroups T^(t) are extended to H by taking T^(t)P*'^ + I - P*^, then these 
extensions are not uniformly exponentially stable and the existence of feedback 
solutions to (R^) on H is not guaranteed. However there is a more subtle 
difficulty regarding verification of the analogue of (2.10) on H (e.g., see 
Theorem 4.3, condition (5.17) of [8]) if the approximate problems are defined 
on H. Even if one has the condition |T*^(t)|j^l^ < Me'“^ with T^(t) extended to H 
as mentioned above, the feedback operators S*^^(t) on H satisfy S^(t)z = z for 
z eH^"^and hence it is not possible to satisfy directly the stability require- 
ment (2.10) on H. In [9 ] this difficulty is handled by taking T*^^(t)z = e ^z 
for z€ But then T(t) and T^(t) are essentially unrelated on H^^. 

In the next two sections we shall see that the version of approximation 
results given in our Theorem 2.2 lends itself to easy verification for certain 
classes of approximation schemes for parabolic systems. 
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3. Convergence of Approximate Riccati Operators for Parabolic Systems 

We use the framework summarized in the previous section to treat the 
optimization of integral quadratic cost functionals for parabolic systems 
of the form given in (1.1). We shall follow the notation introduced in 
section 2 and for our state space H we choose with q a bounded 

domain in possessing a piecewise boundary Unless otherwise 

indicated, all of the function spaces below are to be understood as spaces 
of functions with domain a and range r\ 

For B, D, and Q as given in defining problem (R) of section 2, with 
D > 0 and Q > 0, we consider the regulator problem (R) with the system (2.1) 
for the state y(t) = y(t,*) in H = H^(q) taken as the parabolic system 
n . . n 

y. = I D^(a. .D'^y) + y b.o'y + cy + Bu, t 0, 

^ i,j = l i=l ^ 

(3.1) 


y(0,-) = y(t,-) 


= 0 , 


2 i D 

where u e L (0,«>;U) and D = denotes differentiation with respect to 

<7 X • 

th ^ 

the i spatial variable . 

We make the following standing assumptions on the coefficients in (3.1): 
There exist positive constants y y such that 


I 



1 





for every e G R*^; a.^ = a^^. , and a^.^ G C{u) , b^ G L“(fi), c G L“(j 2 ) 
for every i ,j = 1 , . . . ,n. 

Throughout our discussions the concept of a solution of (3.1) will be that of 
a weak solution (i.e., in the sense of distributional derivatives). 
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We introduce the sesquilinear form 0 : Hq(j 2 ) x Hq(q) C defined by 

0(2, v) = i y a. .D'^zD\ - ( y b.D^z + cz)vldx 

■’a*-i,j=l i = l ^ J 

where c(x) = c(x) - k, with k = k(fj,Y) > 0 determined so that the inequality 

(3.2) Re o(z,z) > [z]^ , z e hJ(j 2) , 

holds for some positive constant c.j independent of z (see[i4, p.l44]). 

Here and throughout I’li and |-| denote the and norms res- 

pectively. Furthermore, to allow use of the theory of sectorial operators 
and sesquilinear forms in discussing the spectra of various operators, we 
assume in defining o that the functions in are complex valued. For 

the sesquilinear form 0 it can be shown (see [14, p. 143] ) that there is also 

a constant C2 = so that 

(3.3) l0(z,v) I 1 C2|z|^ lv|^ 

for all z,v € hJ(j^). Furthermore, it follows from the bounds (3.2), (3.3) 
and well known results on sesquilinear forms (see e.g.,[12, p.lOl]) that 
there exist operators in H { q ) such that 

(3.4) a(z,v) = <-A|^z,v> for z € dom A|^,v e HQ(fi), 

and 

(3.5) o (z , vy = <-A|^v,z> for v € dom A|^, z e Hq(q). 

* 1 

In addition dom Aj^ and dom Aj^ are dense in Hg(fi), 

0**0 

and we have A|^ dom Aj^ = H (fi), A|^ dom Aj^ = H (fi). 
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From (3.2), (3.4) and (3.5), we find that 

2 

Re<A|^z,z> £ -c^|z|.| , z e dom 

(3.6) 

"k O "k 

Re<A|^z,z> 1 -c.||z|^ , z e dom Aj^. 

* 

In view of (3.6) and the range statements above for and Aj^, we may 
invoke standard results from linear semigroup theory [15, p.l6, thm. 4.5] 

k 

to assert that A|^ and A|^ are the infinitesimal generators of linear Cq- 

k 

semigroups Tj^(t) and T|^(t), respectively. (As we have noted above, in 
fact T*(t) = T|^(t)*.) 

We note that the solution semigroup T(t) for (3.1) is given by 

(3.7) T(t) = e''^,^(t) 

with the infinitesimal generator A of T(t) given by A = A|^ + kl and 
dom A = dom A^. Similarly, we have 

(3.8) T(t)* = e'^^T*(t) 

lAr 

with the infinitesimal generator A of T(t) given by A = Aj^ + kl and 
dom A = dom Aj^. We also have [T(t)| £ for t ^ 0. 

Turning next to approximations for (R), we suppose we have a sequence 
of finite dimensional (real) subspaces c HQ(fi), N = 1,2,..., which 
satisfy the approximation condition: 

(Cl): For each z e Hg(i^), there exists an element z^ in H*'* 

such that |z-z^|.| £ f,(N), where e(N) 0 as N -y «. 

N 

We remark that this condition is fulfilled in the event H is chosen 
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as in many classes of finite element approximation schemes [5, Chap. Ill, 3.2J. 

2 

In particular, (Cl) holds for the case where is a rectangle in R and 

N 

the H are the usual linear spans of tensor products of standard one dimen- 
sional piecewise linear splines [16] with mesh size approaching zero as N «>. 

Proceeding in standard fashion, we observe that the restriction of a 

N N N *N 

to H X H defines, in a unique manner, bounded linear operators A|^, 

N 

on H such that 

(3.9) o(z,v) = <-a|^z,v> , z,v e 

and 

N 

(3.10) o(z,v) = <-A|^ v,z> , z,v e H . 

Here A|^^ “ A^ = a|^ + kl, A*^ = a|^ + kl with domains and 

note that A^,A^ generate CQ-semi groups T^^(t), T^(t) on with T^(t) 

the adjoint of T^(t). For the finite dimensional approximating problems 
(R*^) we choose 

(3.11) = P^B, = P^D 

where p''^ H®(jj) is, as in section 2, the canonical orthogonal projec- 
tion. We have thus specified all of the needed entities for the problem (R^) 

of section 2. As we noted previously the trajectories of this problem evolve 
in and hence it is a finite dimensional regulator problem for which com- 
putational techniques are readily available (assuming of course that one has 

N 

made a thoughtful decision in defining the H ). 

We turn to a verification of (H2) of section 2 for the approximations at 
hand. Since it is a trivial matter to see that (Cl) implies that 
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(3.12) -> z, as N -> «>, for z e 

the conditions (H2)-(iii), (H2)-(iv) follow at once from (3.11). We next 

argue that 

(3.13) T|^(t)p '^2 ^ T^(t)z 

(3.14) TjJ(t)Vz T^(t)*z, 

for z e H^(,o), with the convergence uniform in t on bounded subsets of 
[0,«). This taken with (3.7), (3.8) will imply conditions (H2)-(i), (H2)-(ii). 
First we note that from (3.2), (3.9), (3.10) and the fact that e. HQ(fj), 
we have 

(3.15) Re<Aj^z,z> £ -c-| I z| ^ 

(3.16) Re<Aj^ z,z> £-c.||z|^ 

N 

for all z e H . Moreover, we shall demonstrate that the following convergence 
statements hold: 

(3.117) (I-a|J)"''p'^z (I-A^)"''z, z€H°($^), 

(3.18) (I-a[^*)“¥^z -> (I-A*)'^z, z G H°(fi). 

We then may use the Trotter-Kato theorem (see, e.g., the version given in 
[18J) to conclude that (3.15)-(3.18) imply at once the statements (3.13), (3.14). 

Thus we turn to establish (3.17) and (3.18). We shall employ a result 
given in 1 7, p.756. Lemma 3.3], which we state without proof here. 
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Lemma 3.1. There exists a positive constant 6 -^ and a constant 0^ in 
(0 ,ti/ 2) such that 


(3.19) 


+ hli 1 


Xjzp - a(z,z) 


for all z € HQ(fi) and x € C: e.j £ jarg c| £ u|. 

We use this to show that (3.17) holds. For z£ H*^(Q)> define w = 
and = (I~Aj^) ^P^z. Then we have for all z*'* e 




''W,z*^> + a(w,z*^) = <z,z*^>, and 

N N ^ , N Nx N 

<w ,z > + a(w ,z ) = <z,z >. 


N N 

Consequently, defining e = w - w , we find 


N ^ / N N, n 

;e,z > + o(e ,z ) = 0 


for all z^ e Taking X = -1 and z = e*^ in (3.19) - note that 
e^ e h 1($^) - we obtain using this last equation 


N|2 . I N|^ . 

e I + 1 e I ^ £ 6.J 


- e 


a(e^e^ 


= 6-^|-<e ,e +z > - a(e ,e +z )| 

N N N N -N ~N 

for all z e H . Let z = w - w , where w is an approximation for w 
chosen according to (Cl). (Here we again use the fact that w e dom Aj^ c H^Cn) 
We thus obtain the estimate 




< C2e 


(N)«,(|e”| + 


where we have, without loss of generality, assumed that c^ > 1 . This last 



15 


estimate implies e*'* -> 0 in and, in particular, (3.17) holds. 

Turning to (3.18), we recall that 0 ( 2 ,v) = <-A|^v,z> and define 
t: HQ(fi) X C by i(z,v) = a(z,v). Then t satisfies the same 

inequalities (3.2), (3.3) as a. We may therefore verify (3.18) by referring 
to the analysis for (3.17). We summarize our discussions to this point. 

Lemma 3.2 . Let (Cl) hold. Then (H2) obtains with T^(t), and T*^(t) 

defined as in and just above (3.11). 

To use Theorems 2.1, 2.2 of section 2, we shall make the following sta- 
bilizability hypothesis. 

(C2): The pair (A,B) is exponentially stabilizable, i.e., there exists a 

bounded linear operator K: H^(si) ^ U such that the semigroup 
Tg(t) generated by A + BK satisfies |Tg(t)| for 

some positive constants M-j and . 

For a discussion of (C2), we refer the reader to [ 17 J and the references 
given there. 

To make use of the theory of section 2, we need to verify that a certain 
preservation of exponential stabilizability under approximation condition holds 
for our problem. This condition can be stated as: 

(POES): Suppose that condition (C2) holds. Then there exists an integer 

Nq such that for all N ^ Nq the pairs (A^,P*'*B) are uniformly 

exponentially stabilizable by the operator K of (C2), i.e., 

there exist positive constants (independent of N) and 

such that the semigroups T*^(t) generated by A^ + P*^BK satisfy 
N t 

|r^'(t)|<M^e ^ for all N > Nq and t > 0. 
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Before returning to the theoretical results of section 2, we argue that the 
class of approximations for our system (3.1) does indeed satisfy the preserva- 
tion of stabilizability condition (POES). 

Lemma 3.3 . Let (Cl), (C2) hold. Then the approximations defined through 
(3.9), (3.10), (3.11) yield systems that satisfy (POES). 

Proof . We define a sesqui linear form Og: hJ(j^) x Hq(sj) C by og(z,v) = 
a(z,v) + <-BKz,v> + <kz,v> where k is chosen so that 

(3.20) Re Og(z,z) > Cg|z|^ , z e hJ( 5 j) 
and 

(3.21) - ‘=4l^h l''h »z.veHj(fi), 

for some positive constants c^.c^ (recall (3.2), (3.3)). Then arguments 

similar to those in f 7, p.756. Lemmas 3.2, 3.3] can be used to establish 

that the numerical range of ag is contained in a right sector 

S^^ ^ € £ ; |arg(A-Y)| < o| where 0 < 0 < n/2, y real. 

N N 

We next consider the restriction of Og to H x H and, in a manner 

" N " N* 

already discussed, this gives rise to bounded linear operators Ag, Ag on 
such that A^ = A*^ + P^'^BK - k where k = k + k (see the definition of 

D 

~ /.M _ AfJ* 

k in a, k in Og). Indeed og(z,v) = <-Ag z,v>, og(z,v) = <-Ag v,z> 

N '' N* ^ N 

for z,v e H , with Ag the adjoint of Ag. Furthermore, the numerical 

''N/"N*x -f l 

range of A„ (and A^ ) is contained in the left sector S„ = (X: -x G v, 

D B 0 >Y I 9»yJ 

uniformly in N. Thus the numerical range and hence the spectrum (see [11 , 

p.280j) of Ag = A*^ + P^BK = Ag + k are contained in a left sector 

S = S~ + k , uniformly in N. It follows that the set 

0 ,Y 

N 

of all eigenvalues x of Ag with Rex > -6 is bounded, uniformly in N, 
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for any fixed 6. 

Using arguments similar to those behind (3.17), (3.18) - see Lemma 3.1 
and the proof of Lemma 3.2 - it is easily shown that for some r, e r\ ^ > 0 
sufficiently large, we have r. in all the resolvent sets e(Ag), p(A+BK), 

N = 1,2,..., and 

(3.22) (r,-Ag*)"V^z -> (c-(A+BK)*)'^z for z € H^(j^). 

To complete our proof, let us assume that (POES) does not hold even 
though (C2) does. We argue a contradiction. If (POES) does not hold, then 

Ni N 

there exists a sequence N. with N. -> «> and \ an eigenvalue of 

satisfying Re X > -1/j. From our findings on the spectrum of Ap, 

N. 

N = 1,2,..., we know there exists a limit point A of {A with Rex ^0, 

A e S. We shall argue that A is an eigenvalue of A + BK, which is a 
contradiction since (C2) implies Re a £ -w-j for A in the spectrum of A + BK 
(see [ 4, p.32]). 

For convenience, we relabel and drop the subsequent! al notation, assuming 

N '' N N N 

henceforth that A > a, with A an eigenvalue of Ag. Let be an 

eigenvector with |().*^| = 1 and Ag(j)^ = A^())*'* for all N sufficiently large. 

Then we have 




Ab<|) 



and hence 


^-1,,N,N ,Nv _ N 
(Ap-r, ) (a <|) -r,<l> ) - (f) . 


it it 

Let X dom A be arbitrary and put \p = (A+BK-r, ) x- Then 


<A (j) -rj , ((Ag-?) ) P 4<> = <4> ,P ip>. 


(3.23) 
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Using (3.2), one can readily show that the set { 4 )*^} is a bounded set in 

N 

Consequently, there exists a subsequence, again denoted by {<|> }, converging 
strongly in to some nontrivial w e Thus from (3.22), (3.23) 

and (3.12) - a result of (Cl) - we have 

' ( A-r, )w, ( (A+BK-r.) ^) >li> - <w,ijr’ 
or 

<(a-c)w,x> = <w,(A+BK-c) x> 
for all X in dom A . Thus w e dom A and 

<(a-c)w,x> = <{A+BK-c)w,x> 
or 

< (A+BK-A )w»x> - 0 

for all X in dom A . It follows that w is an eigenvector corresponding 

to the eigenvalue a for A + BK and this yields the desired contradiction 

and completes the proof of Lemma 3.3. 

We return finally to a discussion of the convergence theory of section 

2 as it is applied to the specific parabolic system control problem that is 

the focus of the present section. We assume (Cl) and (C2) hold. Then (POES) 

along with Theorem 2.1 yields the existence of nonnegative selfadjoint 

N 

Riccati operators n and n (for N sufficiently large) associated with 
the problems (r) and (r^) in H^(fi) and respecti vely. Since D > 0 
and = P*'^D > 0 on H^, these Riccati operators are unique and furthermore 
(HI) obtains. 

Turning to Theorem 2.2, we first verify that (2.11) holds. 



19 


Faecal 1 that 


(3.24) Ill’ll = sup {<ii^z,z> I z e h’^,|z| = 1} = sup J(z,u^) 

h'^ Izhl 

where is the optimal feedback control (2.8) of (R^). Define, for 

z*^ € with |z^| = 1, the control u^(t) = KT^(t)z^ where T5(t) is 

the semigroup defined in (POES). Then 


1/ M -N, 
J(z ,u ) < 


<D^T^(t)z^T^(t)z^>dt 


:QKT^(t)z^,KT^(t)z^dt 


(|D| + IQIlKl^l M 


2-2 M 

s s 2 


so that from (3.24) we may infer (2.11). To establish (2.10), we first note 

that |S^(t)l for some constants K.j and 3 independent of N. 

H 

This follows from (3.15), (3.16), (2.11) and the fact that S^(t) is generated 
N N -1 * l\i N 

by A - P BQ B P li \ Moreover we have 


0 


-Ds‘^(t)z^,S^(t)z^>dt 


„N N N 

<n z ,z > < 



Since 0 > 0, a theorem of Datko (see [6], [8, p.540, Thm. 2.2]) implies existence 
of positive constants M^ and w, independent of N, such that 


H 

Hence (2.10) of Theorem 2.2 holds. 

Using the convergence results of (2.12), (2.13) it is easy to argue that 
the optimal feedback controls for (R^) converge to that of IR). We summarize 
our findings for the regulator problems for (3.1) in the following theorem. 



Jhcoj'em 3.J. Assume that the subspace approximation condition (Cl) holds 
for H^czHq(c), tiiat the stabi 1 izabi 1 i ty condition (C2) holds for 
(3.1), and that Q ■ 0, D '' 0. Then there exist unique Riccati operators 
and associated with the regulator problems (R) and (R^) on 
and for (3.1) and 


and 


li^P^z >• Hz for z e H^(p.), 
S^(t)P^z ‘ S(t)z for z e H°(^.'); 

u^(t) - u(t). 


with these last two statements holding uniformly in t on compact subsets 
of 10,'"). Here S' (t) and S(t) arc the semigroups generated by 
- P^BQ'^B*p'^n'^ and A - BQ'^B*n., and and u are the 
optimal feedback controls for (r^) and (R), respectively. Moreover, 


-..'t 
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Co n cludinq Remarks 

The conclusions in Theorem 3.1, especially (3.26) and (3.27), are important 

since they reveal that the finite dimensional control laws when employed in 

the systems that we can compute (i.e. , the approximate systems ) allow us to 

anticipate what might happen qualitatively if we used the infinite dimensional 

feedback control s in the original distributed system . However of equal importance 

are findings (which are simple corollaries to the results of section 3) that 

indicate that use of the approximate (readily computed and usually easily 

implemented) controls in the actual distributed system can be expected to produce 

N -1 * N N 

satisfactory performance. More precisely, let U = Q B n P and consider the 
sequence = A - of operators in H®(fi). Then the operators a'^ generate 
semigroups S^(t) which are uniformly exponentially stable and S^(t)z S(t)z, 
uniformly on compact sets in [0,«), zG H^(q), provided, of course, that the 
assumptions of Theorem 3.1 are fulfilled. The uniform exponential stability 
can be established using arguments similar to those in the proof of Lemma 3.3. 

The significance of results for such finite dimensional feedback into the 
original distributed system was noted by Gibson in [9, p. 699]. 

We note that the techniques described in this paper are not restricted to 
parabolic equations of the form (1.1) with distributed control. Indeed as 
can be seen from the arguments in section 3, the essential property required 
for application of these ideas is that the differential equation operator in 
(2.1) be sectorial or, more precisely, that the systems (including feedback) 
generate sesquilinear forms with numerical range in some sector (e.g., see 
the arguments behind Lemmas 3.1, 3.3). Indeed, even though our treatment 
here is concerned with the practically important (in view of the applications 
mentioned in section 1) case of disbributed controls, we recognize that there 
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are important applications where boundary control problems for parabolic equations 
are of primary interest. In some of these applications our treatment and 
techniques are readily used. (We note that the only restriction on B is that 
it be bounded and some boundary control problems are readily transformed to 
the form (2.1)). Furthermore, certain control problems for higher order equations 
can also successfully be treated with the ideas presented in this paper. 

Finally, we note that our approximation approach involves almost no 

N 

restrictions on the subspaces H so that we again can treat a large variety of 

N 

problems. For example, we specifically do not require that H be contained in 
domA (or domAj^) about which we may have only partial information in some 
cases. Thus we may readily employ linear spline approximations with second 
order operators in the framework of our results. Based on our previous efforts 
with spline based approximations in parameter estimation [1, 2] and control 
problems [3J, we are confident that optimism concerning use of splines in the 
present framework is justified. 
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Appendix 


We give here a proof for Theorem 2.2 using in a fundamental way some of 
the results of Gibson. As we have already mentioned, Theorem 2.2 in its present 
form is not given in [s] and in fact his arguments for an analogous result 
appear to contain some technical inaccuracies which we shall attempt to avoid. 

To make our arguments, we need to consider regulator problems on the 
finite intervals [s,t^], - «> < s < t^, with a weighting operator G for the 
final state y(t^). We assume throughout that A generates the Cq semigroup 
T(t) on H, that D, Q, G are selfadjoint with D > 0, Q > 0, G ^ 0, and Bg C(U,H). 
The finite interval problems are given by: 


(R,t^) Minimize J(s,y(s),u) = <Gy(t^-) .y(t|:)> 

+ / {- Dy(t) ,y(t)> + <Qu(t) ,u(t)>}dt 

s 

t 

subject to y(t) = T(t-s)y(s) + / T(t-a)Bu(o)do for s _< t < t^ . 

s 


Under our assumptions a unique nonnegative selfadjoint Riccati operator ri^ 
can be associated with (R,t^). That is, is the unique nonnegative selfadjoint 
solution of the integral Riccati equation for z 6 H 


n. 5 (t)z - T(t^-t)*GT(tf-t)z + / T( 


n-t)*[D- n3(n)BQ“^B*n^(n)]T(ri-t)zdn 


with ng( 5 )€; C(H) for s < 5 < t^, (see [ 8 , Theorems 3.1, 3.2 and equation 
(3.28)]). We then have the following limit results. 

T heorem A. l. Assume that the unique nonnegative selfadjoint solution n of 
the A.R.E. (2.3) exists. Let be the unique Riccati operator function 

associated with the problem (R,0). If lim |S(t)z| = 0 for all z€ H where 

t->- 
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where S(t) is generated by A - BQ ^B*n, then 
(A.1) lim n (s)z = nz for all z H. 

$ ■ ■ > - oo 

If moreover G n and there exist positive constants M and such that 

(A. 2) lS(t)l - , t > 0, 

then 

(A. 3) n _< n^(s) < n + for s < 0. 

Proof . If n is the unique nonnegative, selfadjoint solution of (2.3), then 
by the calculations in [ 8, p. 557-558], it is also the unique solution of the 
first integral Riccati equation of [8] on the infinite interval and corresponds 
to the operator of that paper. Theorem A.l then follows directly from 
Theorem 4. 1 0 of [ 8 ]. 

We note that if in addition to the above hypotheses we have D > 0, then 
(A. 2) is satisfied (see Theorem 4.8 of [8]). 

We next recall an approximation result for the finite horizon regulator 
problem (R,t^) in H. Let (H2) hold with operators as in (s'^) given; in addition 
assume G*^ G C(H^), G^ given. To consider a related finite horizon 

problem in H, we define G^ = G T^ and = D^P^ on H. Consider for 
- 00 < s < t^ and y(s) € H given the problem: 

(n'^.t^) Minimize J^(s,y^(s) ,u) = <GV(tf ) ,y^(t^)> 

+ J DV(t),y^(t) ■ + •Qu(t),u(t)>ldt 
s 

subject to y*^(t) = T^(t-s)P^y(s) + / T*^(t-a)B*^u(a)do for s < t t^ . 

s 
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The problem (k ,t^) is considered as a problem in H even though we note that 
y^(t) e for each t so that D%^(t) = D%''*(t) and G*'*y*^(t^) = G%^(t^). We 
denote the unique nonnegative selfadjoint Riccati operator function associated 
with (R*^,0) by (see Theorem 3.2 of [8]). The following is a consequence 
of Theorem 5.1 of [ 8 ]. 


Theorem A. 2 . Let (H2) hold and assume that > Gz for z€F H. Then for 

s < 0 we have 

u'^ ■> u uniformly on [s,0] , 

y^ y uniformly on [s,0] , 

n^(£)z > iig(rjz for z e: H, uniformly in C on [s,0]. 


-N - -N - 

Here u , u, y , y denote optimal controls and trajectories of the problems 
and (R,0), respectively. 


With these preliminaries, we are now prepared to prove Theorem 2.2. 


Proof of Theorem 2. 2 . 

Denote by and s <0, the Riccati operator functions associated with 
(R,0) and (R^,0) in H where we take G = M 2 I, G^ = M^P^ with the constant in 
inequality (2.11). From Theorem A.l applied to each of the problems (R*^,0) on 
with (2.10) (and hence (A. 2) with M = Mp 3 = w) holding we conclude that 

M 

for s < 0 one has on H 


n^(s) < 11 ^ + M^e^‘^^M 2 


This implies that on H we have for s < 0 and each N 
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(A. 4) II P £ 1 II P + M^e M 2 . 

N N . 

Since IT P IS selfadjoint in H, we conclude from (A. 4) and (A.l) that for 
each t > 0 and z C H, there exists r, = c{z,e) in (-«>,0) such that 


(A.6) |l["(f,)p" - li'V 


.N„B, 


for every N = 1,2,..., 


and 


(A.6) jiiz - 11 ,. ('',)z| • 


Therefore we have 


iV'z - „z| 


iVz - i,"(c)p"z| . |„^(c)p''z - il"(Ozl 


(A. 7) 


+ |iir(^)z - n.(c)zi + liiJOz - Hz| 


,in|z| + |ii^^(c)| |p'^z-z| + ]n^(rjz - n^U)z\ + 


But by Theorem A. 2 and the uniform boundedness principle we have |llp(4)| 
uniformly bounded in N and ii|^{c)z ii ( 4 )z. Finally from (H2)(ii) we have 
P^z -> z and thus (A. 7) implies ii'^P^^z riz for every z€ H. Hence (2.12) is 
establ ished. 

From (H2)(iii) and (2.11) it follows that uniformly 

bounded and moreover B^Q~^B^ ii^P^z > BQ"^B nz for each z €- H. Therefore (2.13) 
follows from use of the variation of parameters representations for 
^(t) = S*^(t)z and y(t) = S(t)z and the Gronwall inequality along with (2.8), 
(2.10) and (H2)(i). Finally (2.14) is a consequence of (2.13) and (2.10). 
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